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Abstract 



Evaluations of the event rates relevant to direct search for dark matter neu- 
tralino are presented for a wide range of neutralino masses and for various 
detector materials of preeminent interest. Differential and total rates are ap- 
propriately weighted over the local neutralino density expected on theoretical 
grounds. 
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1. Introduction 

The most natural way of searching for neutrahno dark matter is provided 
by direct experiments where the effects induced in appropriate detectors 
by neutrahno-nucleus elastic scattering may be measured. The aim of this 
paper is to present the evaluation of the event rates relevant to these direct 
measurements and to discuss the experimental sensitivities required in this 
search for it to be successful. 

Our present analysis refers to a wide range for the neutralino masses, 
20 GeV < < 1 TeV, and examines various detector materials of present 
experimental interest: Ge, Nal, Xe, CaF2. The theoretical framework is 
provided by the Minimal SuSy Standard Model (MSSM), implemented with 
a general GUT assumption, as discussed in the companion paper [1] (this 
reference will hereafter be referred to as paper I). The model depends on a 
number of free parameters which, apart from the masses of a few particles 
which come into play, are chosen here to be M2, tan/?. The definition of 
these parameters is recalled in paper I. 

In the following sections the differential and the integrated rates will 
be presented and discussed also in connection with the properties of the 
neutralino relic density. 



2. Event Rates 

For the nuclear recoil spectrum we use the expression 



where 



Ro,i = NT—(7i <v> (2.2) 

The various quantities in Eq.s (2.1)-(2.2) are defined as follows: Eji 
is the nuclear recoil energy, i.e. En = m^g^f^(l — cos9*)/mN {d* is the 
scattering angle defined in the neutralino-nucleus center of mass frame, thn 
is the nucleus mass, rured is the neutralino-nucleus reduced mass and v is 
the relative velocity), E^"'^ = 2m'^^^v'^/mN, F{Eii) denotes the nuclear form 
factor, a is the neutralino-nucleus cross section. The subscript i refers to 
the two cases of coherent and spin-dependent effective interactions. Nt is 



2 



the number density of the target nuclei and p-^ is the local neutralino matter 
density. Finally, I{Eji) is given by 



In Eq.(2.3) f{v) is the velocity distribution of neutralinos in the Galaxy, 
measured in the Earth rest frame, Vmin(-E'it) is given by Vmin(-E'fl) = 
(mjv-E'it/(2m^g^))^/^. The averages appearing in Eq.(2.2)-(2.3) denote av- 
erages over the velocity distribution in the Earth rest frame. An explicit 
formula for I{Efi) in the case of a maxwellian velocity distribution may be 
found in Ref.[2]. 

The differential rates that will be presented in Sect. 6 will be expressed 
in terms of the electron-equivalent energy E^e rather than in terms of Eji. 
These two variables are simply proportional: Eee = QEn, where Q is defined 
as quenching factor. 

In sect. 6 we will also discuss the total rates, obtained by integrating 
dR/dEee over appropriate ranges of Eee above the experimental energy 
threshold. 

3. Neutralino— Nucleus Elastic Cross Sections 

The total cross sections for neutralino-nucleus elastic scattering have 
been evaluated following the procedure discussed in Ref.s[2-3]. 

Here we only summarize some of the main properties. Neutralino-quark 
scattering is described by the amplitudes with Higgs boson exchanges and 
Z boson exchange in the t-channcl and by the amplitudes with squark ex- 
changes in the s- and u-channels. The neutral Higgs bosons considered here 
are the two CP-even bosons: h, H (of masses nih, mn with mn > Wft) and 
the CP -odd one: A (of mass tjia)- 

The relevant properties for these amplitudes are: 1) Higgs boson ex- 
changes contribute a coherent cross section which is only vanishing when 
there is no zino-higgsino mixtiuc in the neutralino composition[4]; 2) Z bo- 
son exchange provides a spin-dependent cross section which takes contribu- 
tion only from the higgsino components of x; 3) squark exchanges contribute 
a coherent cross section (due to zino-higgsino mixing) as well as a spin- 
dependent cross section (due mainly to the gaugino components of x) [5] • 




(2.3) 



In the evaluations of the cross sections presented here, we have adopted 
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for the masses of the Higgs bosons and of the sfermions the various sets of 
values indicated in paper I. In particular we have relaxed here the assumption, 
used in previous papers [2-3], that rrij > 150 GcV, since this stringent lower 
bound, previously reported by the CDF Collaboration [6], is significantly 
weakened by more refined analyses [7,8]. Thus the following two cases have 
been considered for rrij: 1) = 1.2 m^, when > 45 GeV, = 45 
GcV otherwise, except for the mass of the top scalar partners (the only one 
relevant to radiative corrections) which has been taken m = 3 TeV (hereafter 
this choice for mj will simply be referred to as: rrij — 1.2 m^); 2) = 3 
TeV for all sfermions. The Higgs boson mass rrih (taken here as a free 
parameter) is set equal to values close to its experimental lower bound ~ 50 
GeV [9]. 

Let us recall that, whereas the coherent cross section is proportional 
to the square of the mass number of the nucleus, the non-coherent cross 
section brings a spin-dependence that is commonly represented by a factor 
X^J{J + 1). The values used here for this last quantity are taken from [10]; 
other evaluations may be found in Ref. [11]. 

It is worth noticing that, in the model we are considering here, the 
event rates for neutralino detection are largely dominated by coherent effects 
in most regions of the parameter space. In the small domains where spin- 
dependent effects dominate over the coherent ones the total rates are usually 
too low to allow detection. Then in the present paper we do not consider 
explicitly materials enriched in high-spin isotopes, even if these materials are 
very interesting for a search of hypothetical dark matter particles which in- 
teract with matter via substantial spin-dependent interactions. In particular 
for Ge and Xe we consider only natural isotopic compositions (enrichment 
in some isotope would not appreciably modify our results). 

4. Nuclear Form Factors 

Let us turn now to the £^7j-dependcnce introduced in the nuclear recoil 
spectrum by the form factors. These form factors depend sensitively on the 
nature of the effective interaction involved in the neutralino-nucleus scat- 
tering. For the coherent case, we simply employ the standard parametriza- 
tion[12] 

F(i?^) = ^^^e- W . (4.1) 
where =\ q p= 2mNER is the squared three-momentum transfer, s ~ 
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1 fm is the thickness parameter for the nucleus surface, ro = (r^ — Ss^)"*^/^, 
r = 1.2 A^/^ fm and ji(gro) is the spherical Bessel function of index 1. 

The form factor of Eq.(4.1) introduces substantial suppression in the 
recoil spectrum unless qtq << 1. A noticeable reduction in dR/dEji may 

already occur at threshold En = E^^ = El^/Q when ro ^^tunE^^ is not 
small compared to unity. The actual occurrence of this feature depends 
on a few parameters of the detector material: nuclear radius, quenching 
factor, threshold energy E^. The values of these parameters for the nuclei 
considered in this paper are reported in Table 1 [13]. Also the values of 
F'^{E^^) calculated from Eq.(4.1) are given in this same table. We see that 
the effect of the form factor is very moderate for ^^F,^^Na and ^^Ge. In ^^Ge, 
as compared to ^^F and '^^Na, this occurs since the larger nuclear extension 
is compensated by a smaller value of E^^ and a larger Q. The effect of the 
form factor on the shape of the differential spectrum as a function of £^ee 
depends mainly on how close to the threshold energy El^ the first zero of the 
function jiigro) / (qro) occurs (this zero is located at E^^ ^ 10Q/{m]\[rQ)). 
Comparing the values of El]^ and of E^^ reported in Table 1, we see that only 
for ^^^/ (and marginally for ^"^^Xe) this effect is expected to be noticeable. 

In general for the spin-dependent case there are no analytic expressions 
for the form factors. However, numerical analyses have been performed on 
i3iXe [12], ^^Nb [18] and ^^Ge [19]. The general feature is that these form 
factors have a much milder dependence on En as compared to the coher- 
ent ones, because only a few nucleons participate in the nciitralino-nucleus 
scattering in this case. In our evaluations we use the results of Ref.s [12,19] 
for ^^^Xe and ^^Ge. For ^^^/ we assume that the form factor has the same 
behaviour as in the case of ^^^Xe. For ^^ATa and ^^F we employ the analytic 
expression (4.1), since for these light nuclei the suppression due to nuclear 
size is small anyway. 

Table 1 



Nucleus 


Q 


El^ (KeV) 


F'^iE'^) 


i?Oe(KeV) 




0.09 


8.5 


0.84 


378 




0.23 


4 


0.96 


630 




0.07 


8.5 


0.44 


62 


'^Ge 


0.25 


2 


0.86 


61 


127 J 


0.07 


4 


0.05 


7.4 


'■^'Xe 


0.80 


20 


0.28 


80 
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5. Neutralino Local Density 

As for the value of the local neutralino density, p^, to be used in the 
rate of Eq.(2.2), for any point of the model parameter space we take into 
account the relevant value of the neutralino relic density. When ^^h? is 
larger than a minimal (0/i^)inin consistent with the standard value for the 
local dark matter density pi — 0.3 GeV cm~^, then we simply put p-^ — pi. 
When ^l^y? turns out to be less than (0/i^)min, we take 



Here (0/i^)niin is set equal to 0.05. The effect of this procedure of rescal- 
ing the local dark matter density according to the actual neutralino relic 
abundance is manifest in some characteristic features of our results, as is 
shown in the next section. 



6. Results 

6.1 Germanium. We start the presentation of our results by discussing 
the event rates expected for Ge-detectors. In Fig.l we report the integrated 
rate i?(2 KeV < E^^ < 4 KeV) = f2-^^y dE^edR/ dE^e at the representative 
point: tan/3 = 20, mh — 50 GeV, mj = 1.2 m^^, in the form of a scatter 
plot obtained by varying the M2, p, parameters in the ranges: 20 GeV < 
M2 < 6 TeV and 20 GeV < \p\ < 3 TeV. The range of the neutralino mass 

considered here is 20 GeV < < 1 TeV; however, to simplify the 
presentation of our results, in Fig.l and in the following plots for i?, only the 
range 20 GeV < m^^ < 300 GeV is shown (at higher the plots retain the 
same smooth decrease as m,^ increases). 

To illustrate some of the features of the scatter plot in Fig.l, let us 
recall two main points: i) at fixed m-^ the event rate R depends on the model 
parameters only through the product p-^ai (see Eqs. (2.1)-(2.2)), ii) this 
product is given by 



Px = Pi 



(5.1) 



(n/l2) 



p^iTi = piGi when Vt^h"^ > (Q/i^) 



min 



(6.1) 



and 
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Px^i = ^piCTi oc (f2y/j,^)cri oc ^ when fi^/j,^ < (ri/i^)min 

< <7annV > 

(6.2) 

(f^ann and V being the annihilation cross section and the relative velocity; 
see eq.(3.1) of paper I). 

In Fig.l the wide spread in the values of R at fixed rriy^ is due to the 
very large variations of Px^* ^ scan the M2, /x plane along an isomass 
curve. A detailed analysis of our numerical results shows that the maximal 
value of R, at each m^, usually corresponds to a point in the M2, plane 
where fi^/i^ is substantially below (0/i^)jnin- In other words, the maximal 
expected signal is obtained for neutralino configurations which imply use of 
Eq.(6.2) for p^, with values for the rescaling parameter ^ which may also 
be very small. The nature of the neutralino composition that at fixed rriy 
provides the maximum R depends on a very crucial balance between ai and 
(Tann in the expressiou (6.2). Indicatively, the largest event rates in the plot of 
Fig.l entail the following values for ^: ^ ~ 0.01 — 0.02 for ~ 40 — 60 GeV 
(in this case the gaugino-higgsino mixing is substantial), ^ ~ 0.1 for ~ 
60 — 70 GeV (here the composition is mainly of the higgsino type), ^ ~ 0.05 
for ~ 70 — 90 GeV (large gaugino-higgsino mixing) ; above rriy^ ~ 90 GeV 
the neutralino composition is mainly of the gaugino type with ^ which reaches 
values of order of a few tenths. These properties are depicted in Figs. 2-3. 
Fig. 2 shows fl^h'^ as a function of m^; the diamonds denote the values of the 
relic abundances for those configurations that provide the maximal values 
of R. The compositions of these neutralino configurations are shown in the 
M2-IJ, plane of Fig. 3. 

In view of these results we wish to emphasize that, by disregarding the 
neutralino compositions that give low values of Q^h'^ (say, below Qh'^ ~ 0.05) 
on the base that these compositions are not cosmologically interesting, one 
would miss configurations that potentially generate the most sizable signals 
in the direct searches. 

It is worth noticing that the configurations of maximal gaugino-higgsino 
mixing (defined in paper I) provide large values of R (denoted by circles 
in Fig.l) over the whole m^^ range, even if the rescaling effect is always 
substantial for these compositions; in fact in this case the elastic x-nucleus 
coherent cross sections are very large and compensate for the large values of 
(Tann m Eq.(6.2). 

Because of the properties that we have just discussed, the graph of the 
event rate R (Fig.l) exhibits some of the salient features of the plot of the 
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relic abundance Q-^h^ as a function of (see Fig.2). In particular we 
recognize the pronounced dips at ~ 25 GeV and at ~ 45 GeV, i.e. 
at the values where a"ann has an A{h)-pole and a Z-pole, respectively. 

Finally we wish to stress that in the representative point under discussion 
the maximal expected signal: ~ 6 events/ (Kg day), which occurs at ~ 
40 GeV, differs only by a factor 1.5 from the present experimental sensitivity 
-Rexp — 10 events/ (Kg day) [14] (denoted by a horizontal line in our plot of 
Fig.l). In order to illustrate the discovery potential of the direct neutralino 
search by Ge detectors we show in Fig.4 the regions in the M2, plane that 
could be experimentally investigated in case of improvements in sensitivity 
by a factor of 10 (heavy-dotted regions) or by a factor 100 (light-dotted 
regions). 

In Fig. 5 we give a sample of the differential rates dR/dE^f. corresponding 
to the neutralino compositions that, at a given m^, provide a maximal value 
for the integrated rate R, previously defined. The values m^^ = 40, 80, 150 
GeV are displayed in the figure together with the present experimental points 
[14]. 

In Fig. 6 we give the event rate R at the same values of tan/5 and rrih 
as above, but for the second set of values for the sfermion masses, i.e. i^f = 
3 TeV . We notice that in this case the expected maximal rate is even larger 
than in the previous example for > 70 GeV. This increase is due to 
compositions which are more markedly of the gaugino type (rescaling effect 
is small here, due to the suppression induced by the high values of on the 
annihilation cross section). The relevant explorable regions in the M2 — 
plane are shown in Fig. 7. 

Now let us see how the theoretical rates vary as we change the values for 
tan/3 and for nih- In our evaluations the Higgs-exchange graphs have a sub- 
stantial role; then, increasing rrih has a suppression effect on cross sections, 
and in particular on the elastic one. This is mainly due to the dependence on 
the Higgs masses through the propagators (~ l/m\). tan/? enters in some 
of the coupling constants; it generally occurs that cross sections decrease as 
tan/3 decreases. To illustrate the size of these effects, we show in Fig. 8-9 
the integrated event rates for tan/3 = 8, rrih — 50 GeV for the two sets of 
values for and in Fig.s 10-11 the relevant explorable regions in the M2, ji 
plane. In Fig.s 12-13 R is displayed at the representative point tan/3 = 8, 
rrih = 80 GeV for = 1.2 and = 3 TeV respectively. As compared 
to the previous case tan/3 = 8, rrih = 50 GeV, the present event rate R is 
sizably lower, mainly because of the propagator effect ~ l/w^- The regions 
in the M2,fJ. plane which are experimentally explorable, by increasing the 
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sensitivity appropriately, are shown in Fig.s 14-15. 

We do not report here the results for other representative points. We 
only quote that for instance at tan/? = 2, as is expected, the explorable 
regions are smaller than the previous ones. However, a 2-order-of-magnitude 
improvement in sensitivity would still make feasible an exploration in the 
range 20 GeV < < 100 GeV. 

Then, as far as the Ge detectors are concerned, we can conclude that 
improvements of about 2 orders of magnitude in experimental sensitivities 
would allow investigation of the neutralino dark matter in large regions of 
the model parameter space. It is extremely encouraging that this level of 
sensitivities appears to be within reach in the near future [20] . 

Let us turn now to a brief discussion of the other materials mentioned 
before. The theoretical analysis has been performed along the same lines as in 
the case of Ge. For sake of brevity, here we only report the differential and the 
integrated rates for the two representative points: tan/3 = 20, rrih = 50 GeV, 

= 1.2 and tan/3 = 20, ruh = 50 GeV, = 3 TeV. 

6.2 Sodium Iodide. Theoretical expectations for this diatomic mate- 
rial were already presented and compared with experimental values in Ref . [2] 
for neutralino masses in the range 20 GeV < m-)^ < 80 GeV. Here the analy- 
sis is extended up to m^^ = 1 TeV and to much wider regions of the parameter 
space. The results of the integrated rate i?(4.5 KeV < E'ee < 5 KeV) are 
shown in Fig.s 16-17. 

For the case of tan/3 = 20, nih — 50 GeV, nij = 1.2 we also give in 
Fig. 18 the maximal differential rate at three values of m^. We notice that 
the peculiar behaviour of dR/dE^.^ as a function of is due to the diatomic 
nature of Nal. The response of a Nal detector to the neutralino scattering 
is dominated by Iodine up to E'ee ~ 6 KeV and by Sodium at higher E'ee. 
This circumstance depends on the very different coherent form factor effects 
in the two nuclei as is shown by the values reported in Table I. From the 
threshold up to E^^ ~ 6 KeV in Iodine the large suppression due to F'^[E^^) 
is compensated by a large value of its coherent pointlike cross section {A^ 
effect). At higher values of E'ee the Iodine coherent form factor has a fast 
fall-off, because of the zero at E^e = 7.4 KeV, whereas in the whole range of 
Eqq shown in Fig. 18 the Sodium form factor has a rather flat behaviour. 

As discussed in Ref. [2] an improvement of two orders of magnitude in 
the Nal detectors sensitivity is foreseeable in the next few years. This would 
allow to start the exploration of the physics of dark matter neutralinos with 
these detectors. 
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6.3 Xenon. Much experimental activity is now being developed in dark 

matter detectors which employ Xe. We report in Fig. 19-20 the integrated 
rate R{20 KeV < i?ee ^ 22 KeV) in the representative points: tan/3 = 20, 
nih — 50 GeV, mj = 1.2 m-)^ and tan/3 = 20, ruh — 50 GeV, mj = 3 TeV. 
For the first of these points also the difi"erential rate is given (Fig. 21). We 
notice that for Xe the suppression due to the coherent form factor is sizeable 
(see Table 1). However the use of Xe for detecting dark matter neutralinos 
appear to be very promising in view of the experimental features of these 
detectors. 

6.4 Calcium Fluoride. Finally in Fig.s 22-23 we display the integrated 
rate -R(8.5 KeV < -Eee < 9 KeV) for CaF2 in the usual representative points: 
tan/3 = 20, TJih — 50 GeV, mj = 1.2 m;^ and tan/3 = 20, nih ~ 50 GeV, 

= 3 TeV. In Fig. 24 also the difi^erential rate for the first representative 
point is shown. 

7. Conclusions 

We can conclude that the use of different materials and the achievements 
of improved experimental performances are offering very exciting perspec- 
tives for a direct search of neutralino as a dark matter component. 

The theoretical evaluations presented here aim at setting the level of the 
experimental sensitivities required to start the exploration of the physics of 
dark matter neutralinos. The quantitative relations between the experimen- 
tal sensitivity of a particular detector and its discovery potential are easily 
readable in our plots. 

Some final comments on the main theoretical uncertainties are in order 

here: 

i) major improvements in theoretical predictive power require some experi- 
mental informations about the masses of particles to be discovered at accel- 
erators: Higgs bosons, SuSy particles; 

ii) the evaluation of the elastic coherent cross section suffer from an uncer- 
tainty in the estimate of the Higgs-nucleon strength. In Ref.[3] we have 
analyzed the relevance of this effect for the direct neutralino search, by com- 
paring different estimates for the Higgs-nucleon strength [21,22]. A new 
recent study of this fundamental coupling has shown that the problem needs 
further investigation [23]. In the present paper we have employed for the 
Higgs-nucleon coupling the estimate of Ref.[21]. Use of the estimate in 
Ref.[22] would reduce our rates by a factor of 3 approximately; 

iii) the value of {^h?)m.in to be employed for rescaling the local neutralino 
density is obviously somewhat arbitrary. In the present paper we have con- 
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servatively taken (f2/i^)min = 0.05. In view of the small value of which 

appears to be implied by the estimated amount of dark matter in our galaxy, 
our values for (0/i^)min could quite well be reduced by a factor of 2. This 
reduction in (0/i^)min would increase the maximal rate i? by a factor of 2, 
since, as we have seen above, maximal signals occur in general for neutralino 
compositions where the rescaling effect is large. 

* * * 

This work was supported in part by Research Funds of the Ministero 
deirUniversita e della Ricerca Scientifica e Tecnologica. 
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Table Caption 



Table 1. Quenching factors Q, threshold energies E^^, squared values 
of the form factor (4.1) at the threshold F'^{E^) and the values E^^ where 
the first zero of Eq.(4.1) occurs are listed for the nuclei considered in our 
analysis. 
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Figure Captions 

Figure 1. Scatter plot of the integrated rate R{2 KeV < £^ee < 4 KeV) 
for neutralino-Ge interaction as a function of neutralino mass m^, at the 
representative point: tan/? = 20, nih = 50 GeV and nij = 1.2 m^. The 
parameters M2 and are varied in the ranges 20 GeV < M2 < 6 TeV and 
20 GeV < lul < 3 TeV. The horizontal line denotes the present experimental 
limit -Rexp- Filled-circles represent the R values for neutralino compositions 
where gaugino and higgsino components are maximally mixed (0.45 < af + 
al < 0.55). 

Figure 2. Scatter plot of the neutralino relic abundance O^^/i^ as a 
function of neutralino mass m^. The relic abundance is calculated in the 
same representative point as in Fig.l and also with the same variations over 
M2,//. The horizontal line denotes the value (Sl/i^)min = 0.05. Diamonds 
represent the values of relic abundance which give the maximal value of R 
for a given neutralino mass (they correspond to the maxima of R in 
Fig.l). 

Figure 3. Isomass curves and composition lines for neutralino in the 
M2 — fj, plane for tan/5 = 20. Dashed lines are lines of constant neutralino 
mass {m^ = 30 GeV, 100 GeV, 300 GeV and 1 TeV). Solid lines refer 
to constant gaugino fraction fg in the neutralino composition {fg = af + 
02): fg = 0.99, 0.9, 0.5, 0.1 and 0.01. Squares represent the points in the 
M2 — n parameter space which give the maximal integrated rate R for a given 
neutralino mass m^^ (they correspond to the maxima of R in Fig. 1). 

Figure 4. Explorable regions in the M2 — parameter space for 
neutralino-Ge interaction. Parameters are: 20 GeV < M2 < 6 TeV, 

20 GeV < < 3 TeV, tan/3 = 20, ruh = 50 GeV and mj = 1.2 m^. 
Neutralino masses extend up to = 1 TeV. Heavy-dots denote regions 
where the integrated rate R is in the range 0.1 -Rexp < R < Rexp, being -Rexp 
the present experimental bound. In light-dotted regions 0.01 i?exp < -R < 
0.1 i?exp occurs. Lines denote curves of constant neutralino mass: dashed 
lines stand for = 20 GeV, solid lines for = 100 GeV, dot-dashed lines 
for m^j. = 300 GeV and dotted lines for = 1 TeV. 

Figure 5. Nuclear recoil spectrum dR/dE^e for neutralino-Ge interac- 
tion as a function of the electron-equivalent energy E^q. Lines denote the 
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nuclear recoil spectra which give the maximal integrated rate R for a fixed 
neutralino mass: solid line is calculated for = 40 GeV, dashed line for 
= 80 GeV and dot-dashed line for = 150 GeV. M2 and // parameters 
are varied in the ranges 20 GeV < M2 < 6 TeV and 20 GeV < |//| < 3 TeV. 
The other parameters are: tan/? = 20, nih = 50 GeV and = 1.2 m^^. 
Crosses represent the present experimental limit [14] . 

Figure 6. Same as in Fig. 1, except for = 3 TeV instead of 
my- = 1.2 m^^. 

Figure 7. Same as in Fig. 4, except for m^^ = 3 TeV instead of 
= 1.2 m^. 

Figure 8. Same as in Fig.l; here the representative point is: tan/3 = 8, 
nih — 50 GeV and — 1.2 m-^. 

Figure 9. Same as in Fig.l, here the representative point is: tan/3 = 8, 
nih = 50 GeV and rrij = 3 TeV. 

Figure 10. Same as in Fig. 4, here the representative point is: tan/3 — 8, 
m/i = 50 GeV and = 1.2 m-)^. 

Figure 11. Same as in Fig. 4, here the representative point is: tan/3 = 8, 
ruh = 50 GeV and rrif = 3 TeV. 

Figure 12. Same as in Fig.l, here the representative point is: tan/3 = 8, 
rrih = 80 GeV and = 1.2 m^. 

Figure 13. Same as in Fig.l, here the representative point is: tan/3 = 8, 
ruh = 80 GeV and mf = 3 TeV. 

Figure 14. Same as in Fig. 4, here the representative point is: tan/3 = 8, 
mh — 80 GeV and — 1.2 m^^. 

Figure 15. Same as in Fig. 4, here the representative point is: tan/3 = 8, 
nih = 80 GeV and = 3 TeV. 
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Figure 16. Scatter plot of the integrated rate i?(4.5 KeV < E^e < 

5 KeV) for neutralino-A'^a/ interaction. The representative point is: tan (3 = 
20, rrih = 50 GeV and = 1.2 m^; the parameters M2 and /i are varied 
in the ranges 20 GeV < M2 < 6 TeV and 20 GeV < < 3 TeV. Filled- 
circles represent the R values for neutralino compositions where gaugino and 
higgsino components are maximally mixed (0.45 < af + < 0.55). 

Figure 17. Same as in Fig. 16, except for = 3 TeV instead of 
= 1.2 m^^. 

Figure 18. Nuclear recoil spectrum dR/dE^f. for neutralino-ATa/ in- 
teraction as a function of the electron-equivalent energy i^ee- Lines denote 
the nuclear recoil spectra which give the maximal integrated rate for a fixed 
neutralino mass: solid line is calculated for m-^ — 40 GeV, dashed line for 
m-^ = 80 GeV and dot-dashed line for = 150 GeV. M2 and /j, parameters 
are varied in the ranges 20 GeV < M2 < 6 TeV and 20 GeV < |//| < 3 TeV. 
The other parameters are: tan/3 = 20, rrih = 50 GeV and mj = 1.2 m^. 
Crosses represent the present experimental limit [16]. 

Figure 19. Scatter plot of the integrated rate R{20 KeV < E^e < 
22 KeV) for neutralino~Xe interaction. The values for the model parameters 
and the notations are as in Fig.l6. 

Figure 20. Same as in Fig. 19, except for mj = 3 TeV instead of 
= 1.2 m^. 

Figure 21. Nuclear recoil spectrum dR/dE^.^ for neutralino-Xe inter- 
action as a function of the electron-equivalent energy E^^. The values for the 
model parameters and the notations are as in Fig. 18. Here no experimental 
data are available. 

Figure 22. Scatter plot of the integrated rate i?(8.5 KeV < E'ee < 
9 KeV) for neutralino-CaF2 interaction. The values for the model parame- 
ters and the notations are as in Fig. 16. 

Figure 23. Same as in Fig. 22, except for = 3 TeV instead of 
mi — 1.2 mv 



-17- 



Figure 24. Nuclear recoil spectrum dR/dEee for neutralino-CaF2 in- 
teraction as a function of the electron-equivalent energy E^e- The values for 
the model parameters and the notations are as in Fig. 18. Here no experi- 
mental data are available. 



